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Abstract: Aeolian-fluvial interplay erosion regions are subject to intense soil erosion and are of particular 
concern in loess areas of northwestern China. Understanding the composition, distribution, and transport 
processes of eroded sediments in these regions is of considerable scientific significance for controlling 
soil erosion. In this study, based on laboratory rainfall simulation experiments, we analyzed 
rainfall-induced erosion processes on sand-covered loess slopes (SS) with different sand cover patterns 
(including length and thickness) and uncovered loess slopes (LS) to investigate the influences of sand 
cover on erosion processes of loess slopes in case regions of aeolian-fluvial erosion. The grain-size curves 
of eroded sediments were fitted using the Weibull function. Compositions of eroded sediments under 
different sand cover patterns and rainfall intensities were analyzed to explore sediment transport modes of 
SS. The influences of sand cover amount and pattern on erosion processes of loess slopes were also 
discussed. The results show that sand cover on loess slopes influences the proportion of loess erosion and 
that the compositions of eroded sediments vary between SS and LS. Sand cover on loess slopes 
transforms silt erosion into sand erosion by reducing splash erosion and changing the rainfall-induced 
erosion processes. The percentage of eroded sand from SS in the early stage of runoff and sediment 
generation is always higher than that in the late stage. Sand cover on loess slopes aggravates loess erosion, 
not only by adding sand as additional eroded sediments but also by increasing the amount of eroded loess, 
compared with the loess slopes without sand cover. The influence of sand cover pattern on runoff yield 
and the amount of eroded sediments is larger than that of sand cover amount. Furthermore, given the 
same sand cover pattern, a thicker sand cover could increase sand erosion while a thinner sand cover could 
aggravate loess erosion. This difference explains the existence of intense erosion on slopes that are thinly 
covered with sand in regions where aeolian erosion and fluvial erosion interact. 
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1 Introduction 


Erosion is an important factor that influences environmental quality and ecological stability. 
Co-erosion is a complex mechanism caused by various types of erosion, such as aeolian erosion 
and fluvial erosion (Ta et al., 2014; Liu et al., 2016). One typical area containing these two 
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erosion types (i.e., aeolian erosion and fluvial erosion) is the ten small tributaries to the Inner 
Mongolia Reach of the Yellow River where the Hobq Deserts is located (Li et al., 2015). This area 
is severely subjected to aeolian erosion and fluvial erosion (occurring alternately throughout the 
year), and it is the main source of coarse sediments in the Yellow River (Xu, 2014). In the ten 
small tributaries to the Inner Mongolia Reach of the Yellow River, the yield of eroded sediments 
and the length of erosion time are both increased by the superimposed effect of erosive power and 
erosion time, resulting in greater environmental and ecological influences. During the fluvial 
erosion period, the short-duration, high-intensity erosion prevails, and dryland fluvial systems 
often lead to intense soil erosion and sedimentation (Kocurek, 1998; Banham and Mountney, 
2014). Therefore, understanding the mechanisms of aeolian erosion and fluvial erosion is of great 
importance to control and prevent the loss of soils. 

In the ten small tributaries to the Inner Mongolia Reach of the Yellow River, the sand-covered 
loess slopes (SS) are widely distributed in the upstream areas (Zhang et al., 2015). The thickness 
of sand layer is ~1.5 cm on average, and the maximum slope gradient of the SS is confirmed as 
~12° in gully heads from both field investigations and ArcGIS analysis using a 90-m resolution 
digital elevation model. Zhang et al. (2015) found that the upstream watershed section is the main 
source of eroded sediments and that the erosion dynamics are controlled by the fluvial processes. 
In later work, Yao et al. (2016) pointed out that the transported sediments in the Xiliugou Basin of 
the upper Yellow River (i.e., a region where aeolian erosion and fluvial erosion interact) is 
generated from the SS; they, therefore, concluded that SS are a key source of eroded sediments. 
This means that research on SS is important to understand the mechanisms of aeolian-fluvial 
erosion processes in the ten small tributaries to the Inner Mongolia Reach of the Yellow River. 

Laboratory rainfall simulation experiments are one crucial approach for the study of the 
mechanisms of erosion on slopes. As previous studies have shown, during short-duration rainfall 
events, the rainfall intensity has a strong influence on runoff and sediment generation on slopes 
(Huang et al., 2016). Moreover, variations in the underlying surface conditions, such as plant 
cover (Garcia-Estringana et al., 2010), and sand, stone or gravel cover (Abrahams et al., 2001; 
Mandal et al., 2005; Turowski et al., 2010; Zhang et al., 2015), can also influence runoff 
generation and total eroded sediments by changing soil erosion processes (Wang and Shi, 2015). 
Compared with uncovered slopes, gravel cover on slopes can alter the characteristics of eroded 
sediments (Turowski et al., 2010). Stones resting on the soil surface can cause a pronounced 
reduction of runoff and sediment erosion and an increase of infiltration compared with those 
buried in the surface layer (Sun et al., 2002; Mandal et al., 2005). The sand-covered slopes can 
generate higher amount of eroded sediments and greater median grain size of eroded sediments 
than the uncovered slopes (Xu et al., 2015). To date, few studies have systematically examined 
the erosion processes of SS at different rainfall intensities in aeolian-fluvial erosion regions. 

The objective of this research was to investigate the influences of sand cover on erosion 
processes of loess slopes using rainfall simulation experiments. The composition, distribution, 
and transport processes of eroded sediments from SS and uncovered loess slopes (LS) were 
explored. Specifically, we fitted the grain-size curves of eroded sediments, and obtained the 
composition and source of eroded sediments from SS and LS based on grain-size distribution 
curves. We also described the variability characteristics of eroded sediments during runoff and 
sediment generation, and analyzed the influences of rainfall intensity, sand cover amount, and 
sand cover pattern on erosion processes. It is our hope that this research will help to elucidate the 
influences of sand cover on erosion processes of loess slopes and to provide a rationale for further 
in-depth studies of the co-erosion processes in the aeolian-fluvial interplay erosion regions. 


2 Materials and methods 


2.1 Study area 


The loess samples used in the rainfall simulation experiments were collected from the upper 
watershed of Dongliu Gully, which is one of the ten small tributaries to the Inner Mongolia Reach 
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of the Yellow River. The aeolian sand samples were collected from the Hobq Deserts, which is 
crossing the middle watershed of the Dongliu Gully. The Dongliu Gully watershed has a typical 
continental climate, with long, cold winters and short, hot summers. Runoff is primarily recharged 
by floods following heavy rain, with rainfall occurring mainly in the form of storms. Fluvial 
erosion in the study area is dominated by rainstorm-induced floods. Grain-size characteristics of 
the loess and sand in the study area are shown in Table 1. 


Table 1 Grain-size characteristics of the experimental loess and sand 


Type dso (mm) Clay (%) Silt (%) Sand (%) 
Loess 0.010 13.20 82.53 4.27 
Sand 0.418 0.09 3.30 96.61 


Note: dso is the medium value of the grain-size distribution. 


2.2 Experimental design and set-up 


Before the rainfall simulation experiments, 9 loess samples and 9 sand samples were collected for 
the analyses of the basic sediment grain size. Grain-size analysis of the loess and sand samples 
was conducted using a Malvern Mastersizer 2000 laser grain-size analyzer (Malvern, UK). The 
sediment grain-size curves of these samples were averaged and taken as the background 
grain-size curves. 

All collected loess and sand samples used to rainfall simulation experiments were sieved with a 
10-mm sieve before they were packed into the soil bin. Rainfall simulation experiments were 
conducted at the water and soil erosion laboratory of Xi’an University of Technology during 
July-September in 2013. We used a slope test apparatus to perform laboratory simulation 
experiments. The experimental soil bin comprised a mobile steel container with 13.0 mx1.0 
mx0.6 m dimensions. The experimental slope gradient was set at 12°. The reasons are as follows. 
First, it has been proven that the maximum slope gradient of the SS is ~12° in the gully head of 
the ten small tributaries to the Inner Mongolia Reach of the Yellow River (Zhang et al., 2015). 
Second, Foster et al. (1984) concluded that the critical slope gradient of erosion is less than 25°. 
Third, Zheng et al. (2005) and Shen et al. (2016) found that the erosion characteristics at slope 
gradients of 3°-12° are different from those at slope gradients of 12°-25°. We therefore 
considered 12° to be the critical slope gradient in this study and used it as our experimental slope 
gradient. 

The schematic diagram and stereogram of the simulated rainfall system used in this study are 
shown in Figure 1. Specifically, the lowermost layer of the soil bin was filled with 0.1-m-thick 
coarse gravel wrapped in gauze to ensure water permeability, and a 0.2-m-thick coarse gravel was 
laid at the bottom of the bin to keep the soil drainage conditions of the test close to those of a 
natural slope. The middle layer of the soil bin comprised a 0.25-m-thick experimental loess soil 
samples, with mean bulk density of 1.3 g/cm? and water content of 20%. The uppermost layer of 
the soil bin encompassed sand cover with different patterns (including thickness and length). The 
uppermost layer of the soil bin without sand cover (i.e., LS) was considered as control. According 
to the results of Ta et al. (2014), we designed three sand cover thicknesses (0.5, 1.0, and 1.5 cm) 
that are relatively common in aeolian-fluvial erosion regions. The total length of loess slope is 
13.0 m. The lengths of sand cover on loess slopes were designed as 13.0, 9.9, 6.6, and 3.3 m, 
corresponding to the whole slope, about three-quarters of slope, about one-half slope, and about a 
quarter of slope, respectively. It should be noted that if the sand cover length in any experimental 
test was less than 13.0 m, then this length was calculated upwards from the bottom of the soil bin. 
Once the sand cover was established within the test apparatus, two points were chosen randomly 
within each meter to ensure that the sand layer thickness was constant. Therefore, a total of 26 
randomly chosen points were selected for this study. 

A side sprinkler device was used to simulate rainfall from an effective height of 4 m. We 
applied three levels of rainfall intensity (1.0, 1.5, and 2.0 mm/min) based on the rainfall 
characteristics of the study area. For rainfall intensity levels of 1.0 and 1.5 mm/min, the rainfall 
duration was set to 30 min; while for rainfall intensity level of 2 mm/min, the rainfall duration 
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was set to 15 min because of a water feeding issue. Finally, a total of 11 experimental tests were 
performed (Table 2), with three replicates for each. The results presented in the following are the 
average of the three replicates. 
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Fig. 1 (a) Schematic diagram and (b) stereogram of the simulated rainfall system 


Table 2 Design of experimental tests for the rainfall simulation experiments 


Test No. Sand cover length (m) Sand cover thickness (cm) Rainfall intensity (mm/min) 
LS-1.5 0.0 0.0 1.5 
LS-1.0 0.0 0.0 1.0 
LS-2.0 0.0 0.0 2.0 

SS-13.0-1.0-2.0 13.0 1.0 2.0 
SS-13.0-1.0-1.5 13.0 1.0 1.5 
SS-13.0-1.0-1.0 13.0 1.0 1.0 
SS-13.0-0.5-1.5 13.0 0.5 1.5 
SS-13.0-1.5-1.5 13.0 1.5 1.5 
SS-3.3-1.0-1.5 3.3 1.0 1.5 
SS-9.9-1.0-1.5 9.9 1.0 1.5 
SS-6.6-1.0-1.5 6.6 1.0 1.5 


Note: LS, uncovered loess slopes; SS, sand-covered loess slopes. The test number was expressed as loess slope type-sand cover length 


(m)-sand cover thickness (cm)-rainfall intensity (mm/min). 


2.3 Experimental process and data analysis 


The water was sprinkled evenly across the slopes to ensure that water infiltration was closer to the 
natural soil water distribution, then the simulated slopes were allowed to stand for 24 h before the 
simulated rainfall events. The water content in the soil was controlled at 20%. Before the 
experiments, rainfall intensity was calibrated to test its uniformity. During the experiments, the 
amount of runoff and eroded sediments per minute was sampled using a pycnometer to calculate 
the sediment concentration. The remaining runoff and eroded sediments were collected using a 
bucket. The water and eroded sediments in the pycnometer were added to the bucket after the 
calculation of sediment concentration. Finally, the water and sediments in the bucket were 
weighed to calculate the amount of water and of eroded sediments produced in 1 min. The amount 
of eroded sediments per minute in the runoff was measured using an oven-drying method. All 
eroded sediments were precipitated and placed in mess tins in a drying oven at 105°C for 8 h (Shi 
et al., 2012). On both sides of the soil bin, 0.5-m-high splash erosion plates were set upright to 
collect the splash-eroded sediments. All samples were sieved with a 2-mm sieve. Then a Malvern 
Mastersizer 2000 laser grain-size analyzer (Malvern Instruments, Ltd., UK) was used to 
determine the grain size of eroded sediments. 
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Grain characteristics, including size and distribution, are related to transport media, transport 
modes, and dynamic erosion conditions. The combination of eroded sediments from multiple 
sources leads to variations in the number and position of peaks in the grain-size distribution curve, 
and therefore, variations in the curve can reflect the changes in sediment composition (Weltje and 
Prins, 2007). If a grain-size distribution curve follows a polymodal distribution, then it can be 
concluded that the distribution is made up of multiple, superposed standard unimodal functions. 
In other words, the polymodal grain-size distributions of loess and sand sediments are resulted 
from multiple single-fraction. 

Three distribution functions are fitted to the sediment grain-size distribution: normal (Sun et al., 
2001), lognormal (Xiao et al., 2013), and Weibull (Sun et al., 2002; Zhang et al., 2015). However, 
as the normal function has a zero skewness coefficient, it is not a suitable fit for the variety of 
grain-size distribution (Sun et al., 2001). Sun et al. (2002) noted that the Weibull function fits 
better than its lognormal counterpart in the case of aeolian and hydraulic sediments; Zhang et al. 
(2015) also showed that the Weibull function encompasses an appropriate distribution when 
characterizing sediment from aeolian and fluvial erosion, providing satisfactory results for 
sand-loess mixtures. Therefore, we selected the Weibull function to describe the sediment 
grain-size distribution in this study. The equation of the Weibull function is as follows (Sun et al., 
2002; Zhang et al., 2015): 


k-1 
f(x4,k)= ka) e764} Gf x20), (1) 
f (x;4,k) =0 (if x<0), (2) 


where x is a random variable that represents the grain-size fraction of the sediment; 4 is a scale 
parameter; and k is a shape parameter. Each unimodal component can be described with sufficient 
accuracy by the Weibull function, and the distributions of polymodal sediments that consist of 
several unimodal fractions are considered to be a mixed grain-size distribution of the unimodal 
fractions. Assuming that the distribution function f(x; A, k) of polymodal sediments has four 
possible sources, i.e., fi(x1; 21, k1), fo(x2; A2, k2), fa(x3; A3, k3) and fa(xa; A4, k4), then the following 
relationship is satisfied: 

f(t A,k) = af (x15 Ay ky) + Bf (X25 Ap sky) + Of (535 ks) + dfa (X43 Ags ky) 

a+b+c+d=1 


The Weibull distribution parameters of the unimodal fractions were first estimated by 
minimizing the residual differences between the measured and fitted grain-size curve of the basic 
sediment samples. The percentage parameters, i.e., a, b, c and d, for the polymodal sediments were 
then determined by minimizing the residual differences between the measured and fitted grain-size 
curve of the eroded sediments from LS and SS. The calculation is conducted by R compiler. 


i (3) 


3 Results 


3.1 Basic characteristics of sediment grain size 


As shown in Figure 2, the basic grain-size distribution curves of the loess and sand used in the 
rainfall simulation experiments exhibited different bimodal distributions. Using the two grain-size 
curves of the basic loess and sand samples as objective functions, we applied two Weibull 
functions for fitting without constraints, with the minimum fitting residuals as the fitting results of 
the two functions. In this case, C1 refers to the fine loess grains (accounting for 36.2% of the 
loess grain), which are mainly composed of clay and have a dso of 0.003 mm. Similarly, C2 refers 
to the coarse loess grains (accounting for 63.8% of the loess grain), which are mainly composed 
of fine silt particles and have a dso of 0.016 mm. C3 refers to the fine sand grains (accounting for 
2.8% of the sand grain), which are mainly composed of coarse silt particles and have a dso of 
0.041 mm. Finally, C4 refers to the coarse sand grains (accounting for 97.2% of the sand grain), 
which are mainly composed of medium sand particles and have a dso of 0.373 mm. 
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Fig. 2 Model fitting results for the basic grain-size distributions of loess (a) and sand (b) used in the rainfall 
simulation experiments. C1, the fine loess grains; C2, the coarse loess grains; C3, the fine sand grains; C4, the 
coarse sand grains. 


The grain-size characteristics of the four fractions (C1, C2, C3 and C4) are shown in Table 3. 
The dso values of the four fractions increased in turn. The sorting coefficients were different 
between loess (C1 and C2) and sand (C3 and C4), with the sorting behavior of sand being better 
than that of loess. All four distribution curves showed positive skewness, and their kurtoses were 
moderate with little differences. To sum up, the dso, kurtosis, skewness, and sorting coefficient 
values varied with the fraction (C1—C4), reflecting that their differences were quite large and the 
large differences were the preconditions for the gradation curve (grain-size distribution) analyses. 


Table 3 Grain-size characteristics of the four fractions 


Fraction dso (mm) Sorting coefficient Skewness Kurtosis 
Cl 0.003 3.83 0.70 1.00 
C2 0.016 3.07 0.10 1.02 
C3 0.041 1.92 0.12 1.04 
C4 0.373 1.85 0.13 1.04 


Note: C1, the fine loess grains; C2, the coarse loess grains; C3, the fine sand grains; C4, the coarse sand grains. 


3.2 Composition of eroded sediments from SS and LS 


Using the grain-size distributions of eroded sediments in runoff per minute during the 11 
experimental tests as the objective functions, we fitted the four fractions (C1—C4) using least 
squares to calculate the sand and loess compositions in the total eroded sediments from SS and LS 
under different simulated conditions. From this analysis, C1+C2 fraction was attributed to loess 
and C3+C4 fraction was attributed to sand. The mean compositional percentages of eroded 
sediments in all experimental tests are shown in Table 4. 

As demonstrated in Table 4, regardless of the loess slopes with or without sand cover, the main 
eroded sediments were loess, with the percentage of C1+C2 fraction accounting for >83% of the 
total eroded sediments on average. Furthermore, C2 fraction was the predominant and 
respectively accounted for 86% and 68% of the total eroded sediments from LS and SS, 
indicating that variations in the underlying surface conditions could influence the composition of 
eroded sediments. The sand percentage in the total eroded sediments varied with the sand cover 
patterns and rainfall intensities, with the percentage of C3+C4 fraction varying over a range from 
9% to 25%. Furthermore, C4 fraction was the predominant of eroded sand. Loess slopes with 
sand cover could influence the composition of eroded loess. Specifically, the mean ratio of C1 
fraction to C2 fraction for LS was 0.162, while it was 0.221 for SS. Hence, sand cover generally 
reduced the proportion of C2 erosion on loess slopes. 

In this study, the splash-eroded sediments from LS was collected using splash erosion plates 
and then further analyzed. The results showed that the dso of the splash-eroded sediments was 
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0.015 mm. Grain-size distribution curves of the splash-eroded sediments and C2 fraction from LS 
are shown in Figure 3. Generally speaking, the two distribution curves were very similar in shape, 
but the sorting coefficient of the splash-eroded sediments (3.56) was higher than that of C2 
fraction (3.07), indicating that the splash-eroded sediments were generated through further sorting 
of C2 from LS. 


Table 4 Mean compositional percentage of eroded sediments in different experimental tests 


Compositional percentage (%) 


Test No. 
Cl fraction C2 fraction C3 fraction C4 fraction 

LS-1.5 15 85 0 0 
LS-1.0 15 85 0 0 
LS-2.0 13 87 0 0 
Mean 14 86 0 0 
SS-13.0-1.0-2.0 13 71 0 16 
SS-13.0-1.0-1.5 11 67 0 22 
SS-13.0-1.0-1.0 18 65 2 15 
SS-13.0-0.5-1.5 14 76 1 9 
SS-13.0-1.5-1.5 15 63 0 22 
SS-3.3-1.0-1.5 15 68 1 16 
SS-9.9-1.0-1.5 21 70 0 9 
SS-6.6-1.0-1.5 15 60 1 24 
Mean 15 68 1 17 
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Fig. 3 Comparison of grain-size curves between splash-eroded sediments and C2 fraction from LS (uncovered 
loess slopes) 


The intensity of splash erosion is related to the clay and silt percentage of soil (Saedi et al., 
2016) and the soil aggregate properties (Abu-Hamdeh et al., 2006). In this study, the total clay 
and silt content in the sand layer was only 3.39%; in contrast, the clay and silt content in the loess 
layer was 95.73% (Table 1). The lower clay and silt content in the sand layer relative to the 
underlying surface conditions led to the extremely low splash erosion intensity of SS in the initial 
of runoff and sediment generation. This suggests that sand cover on the loess slopes could change 
the rainfall-induced slope erosion process by reducing splash erosion. 


3.3 Variability characteristics of eroded sediments during runoff and sediment generation 


Although all experimental results demonstrated that eroded sediments were dominated by C2 
fraction throughout the runoff and sediment generation, the compositional percentages of eroded 
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sediments changed to various extents. The peak occurrence time, peak sand percentage and mean 
sand percentage of eroded sand during the first 5 min of runoff and sediment generation in 
different experimental tests are shown in Table 5. One noteworthy result is that it is possible to 
define three stages to describe the process of runoff and sediment generation. The first (early) 
stage occurred from the 1* to the 5" min of runoff and sediment generation at rainfall intensities 
of 1.0 and 1.5 mm/min. The second (intermediate) stage encompassed the period from the 6" to 
the 20" min of runoff and sediment generation at rainfall intensities of 1.0 and 1.5 mm/min, and 
from the 6" to the 10" min of runoff and sediment generation at the rainfall intensity of 2 mm/min. 
The third (late) stage occurred from the 21" to the 30" min of runoff and sediment generation at 
rainfall intensities of 1.0 and 1.5 mm/min. It should be noted that the late stage at the rainfall 
intensity of 2 mm/min occurred from the 11" to the 15" min of runoff and sediment generation. 


Table 5 Peak occurrence time, peak sand percentage, and mean sand percentage of eroded sand during the first 
5 min of runoff and sediment generation in different experimental tests 


Test No. Peak occurrence Peak percentage of C3+C4 Mean percentage of C3+C4 fraction in the first 5 
time (min) fraction (%) min (%) 
LS-1.5 0 0 0 
LS-1.0 0 0 0 
LS-2.0 0 0 0 
SS-13.0-1.0-2.0 3 55 34 
SS-13.0-1.0-1.5 3 91 69 
SS-13.0-1.0-1.0 5 69 46 
SS-13.0-0.5-1.5 2 65 43 
SS-13.0-1.5-1.5 3 94 60 
SS-3.3-1.0-1.5 5,13 39 19 
SS-9.9-1.0-1.5 3 52 19 
SS-6.6-1.0-1.5 4 73 34 


Note: Experimental test SS-3.3-1.0-1.5 produced two peaks of similar magnitude, occurring at the 5 and 13" min of runoff and 
sediment generation. 


The results presented in Table 5 show that the peak sand (C3+C4 fraction) percentage in eroded 
sand was generally observed in the period of the 314-5 min, with the percentage value ranging 
from 39% to 94%. Just one experiment (SS-3.3-1.0-1.5) produced two peaks of similar magnitude, 
occurring at the 5“ and 13" min. This might because the amount and length of sand cover in this 
experiment were the smallest among the rainfall experiments, resulting in small and unstable 
peaks. The mean percentage of eroded sand from SS within the first 5 min of runoff and sediment 
generation ranged from 19% to 69% (Table 5), while the mean sand percentage of eroded 
sediments from SS in the whole process was 17% (Table 4). These results show that sand erosion 
mainly occurs during the early stage of runoff and sediment generation. 

To study the influence of rainfall intensity on the composition of eroded sediments, we 
compared the percentage variations of the main fraction (C2) in eroded sediments from SS and LS 
along with erosion time under three rainfall intensity levels (1.0, 1.5, 2.0 mm/min), and the results 
are presented in Figure 4 and Table 6. For LS, the percentage of C2 in eroded sediments under the 
three rainfall intensity levels ranged from 81% to 89% (Fig. 4), and all coefficient of variation 
(Cy) values for C2 over the whole rainfall event were relatively small (Table 6). These results 
indicate that the composition of eroded sediments from LS did not change significantly with 
rainfall duration and rainfall intensity. However, the composition of eroded sediments from SS 
changed markedly with runoff and sediment generation duration and rainfall intensity. The Cy 
values for C2 were the largest, and C2 changed most notably at the rainfall intensity of 1.5 
mm/min, whereas the composition of eroded sediments changed only slightly at the highest 
rainfall intensity (2 mm/min). 

As shown in Figure 4, the percentage of C2 in eroded sediments from SS during the early stage 
of runoff and sediment generation varied dramatically at all rainfall intensity levels; however, it 
changed slightly during the late stage. By comparison, the percentage of C2 in eroded sediments 
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from SS during the whole runoff and sediment generation changed slightly. These results mean 
that sand cover had influence on the erosion processes of loess slopes, and the changes mainly 
occurred in the early stage of runoff and sediment generation. As rainfall proceeded, stable 
erosion gullies formed, and the characteristics of eroded sediments from SS gradually became 
more similar to those from LS. In other words, the compositions of eroded sediments from SS and 
LS became more similar. 
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Fig. 4 Percentage variations of C2 in eroded sediments from SS (sand-covered loess slopes) and LS (uncovered 
loess slopes) along with erosion time under rainfall intensities of (a) 1.0 mm/min, (b) 1.5 mm/min, and (c) 2.0 
mm/min 


Table 6 Coefficient of variation (Cv) values for C2 in eroded sediments from SS and LS during different stages 
of runoff and sediment generation in the selected experimental tests 


Runoff and sediment generation LS-1.0 SS-13.0-1.0-1.0 LS-1.5 SS-13.0-1.0-1.5 LS-2.0 SS-13.0-1.0-2.0 


Whole duration 0.02 0.26 0.04 0.33 0.01 0.23 
Early stage 0.01 0.43 0.05 0.95 0.01 0.37 
Intermediate stage 0.02 0.22 0.04 0.11 0.01 0.05 
Late stage 0.02 0.04 0.04 0.09 0.01 0.03 


4 Discussion 


4.1 Influences of sand cover amount on erosion processes of loess slopes 


Changes in ambient conditions can influence the erosion processes of slopes (Huang et al., 2016). 
The influences of rainfall intensity and sand cover pattern on erosion processes of slopes can be 
reflected by the changes in runoff yield and amount of eroded sediments (Rodriguez-Lopez et al., 
2012). In this study, the amount of sand cover and corresponding total amount of eroded sediments 
in different experimental tests are shown in Table 7. These data could be used to further explore the 
influences of rainfall intensity and sand cover pattern on erosion processes of loess slopes. 
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Table 7 Amounts of sand cover and total eroded sediments recorded in different experimental tests 


Test No. Amount Fi hdieover Total amount of eroded sediments (kg) ee eta 
LS-1.0 0.00 11.91 
SS-13.0-1.0-1.0 344.50 275.90 263.99 
LS-1.5 0.00 38.69 
SS-13.0-1.0-1.5 344.50 513.39 474.70 
SS-13.0-0.5-1.5 172.25 585.83 547.14 
SS-13.0-1.5-1.5 516.75 855.76 817.07 
SS-3.3-1.0-1.5 87.45 94.79 56.10 
SS-9.9-1.0-1.5 262.35 285.71 247.02 
SS-6.6-1.0-1.5 172.25 313.10 274.41 
LS-2.0 0.00 214.09 
SS-13.0-1.0-2.0 344.50 540.83 362.74 


Note: The increment of eroded sediments from SS relative to from LS was calculated by the difference between the total amount of 
eroded sediments from SS and the total amount of eroded sediments from LS at the same rainfall intensity. 


Given the same rainfall intensity, the amounts of eroded sediments from SS were larger than 
those from LS. The increments of eroded sediments for SS-13.0-1.0-1.5, SS-13.0-0.5-1.5, 
SS-13.0-1.5-1.5, SS-6.6-1.0-1.5, and SS-13.0-1.0-2.0 were all larger than the corresponding 
amounts of sand cover, demonstrating that sand cover not only provided additional eroded 
sediments for loess slopes but also exacerbated the erosion of loess slopes to a certain extent. 
There are also some other studies demonstrated that the sand layer on the loess slopes could 
significantly influence the processes and mechanisms of runoff and sediment generation (Xu et al., 
2015; Zhang et al., 2017a). The sand cover on loess slopes may destroy the stability of the water 
flow, further aggravating soil erosion of loess slopes (Zhang et al., 2017b). 

To analyze the increments of eroded sediments from SS relative to from LS, we calculated the 
amount of sand cover in eroded sediments per minute by multiplied each sand cover proportion in 
eroded sediments by the corresponding amount of eroded sediments per minute. Variations in the 
amounts of C1+C2 and C3+C4 in the early and late stages of runoff and sediment generation 
under different amounts of sand cover are shown in Figure 5. The amount of C1+C2 represented 
the mass of eroded loess, and the amount of C3+C4 represented the mass of eroded sand. 


60 


180 


160 f -*-C1+C2 in the early stage (b) 
E) —— C1+C2 in the late stage B 504 + C3+C4 in the early 
z 140 \ = stage 
8 120 @ 40 | —— C3+C4 in the late 
an stage 
3 100 3 
Ò S 304 
5 80 8 
S © 
ž 60 E 20 4 
£ 40 E 
10 4 
20 
0# f 0 7 : z - - 
0 100 200 300 400 500 600 0 100 200 300 400 500 600 
Amount of sand cover (kg) Amount of sand cover (kg) 


Fig. 5 Variations in amounts of eroded loess (a) and eroded sand (b) in the early and late stages of runoff and 
sediment generation under different amounts of sand cover 


From Figure 5 it can be concluded that the amount of eroded loess (C1+C2) in the early stage 
of runoff and sediment generation was always smaller than that in the late stage, especially for the 
172.25 kg amount of sand cover. By contrast, the amount of eroded sand (C3+C4) in the early 
stage was always larger than that in the late stage, with an exception of the 172.25 kg amount of 
sand cover. In the early stage, the amount of eroded sand (C3+C4) in the SS-13.0-0.5-1.5 case 
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(172.25 kg amount of sand cover) was larger than that in the SS-6.6-1.0-1.5 case (172.25 kg 
amount of sand cover); while in the late stage, the amount of eroded sand (C3+C4) was larger in 
the SS-6.6-1.0-1.5 case than in the SS-13.0-0.5-1.5 case. These results indicate that sand cover 
pattern has considerable influence on erosion processes of loess slopes. 

4.2 Influences of sand cover pattern on erosion processes of loess slopes 

In this study, we conducted paired comparison analyses for experimental tests of SS-13.0-1.0-1.5, 
SS-13.0-0.5-1.5, SS-6.6-1.0-1.5, and LS-1.5 under the rainfall intensity of 1.5 mm/min to study 
the influences of sand cover pattern on erosion processes of loess slopes (Table 8). It should be 
noted that tests SS-13.0-1.0-1.5 and SS-13.0-0.5-1.5 had the same sand cover length but different 
sand cover amounts and thicknesses, while tests SS-13.0-0.5-1.5 and SS-6.6-1.0-1.5 had the same 
amount of sand cover but different sand cover patterns (thickness and length). 


Table 8 Amounts of sand cover, total eroded sediments, eroded sand and eroded loess as well as runoff yield in 
the selected four experimental tests 


Test No. Amount of sand Amount of total eroded Amount of total Amount of total Runoff yield 
cover (kg) sediments (kg) eroded sand (kg) eroded loess (kg) (m?) 
SS-13.0-1.0-1.5 344.50 513.39 110.17 403.22 455.47 
SS-13.0-0.5-1.5 172.25 585.83 44.69 541.14 501.99 
SS-6.6-1.0-1.5 172.25 313.10 68.93 244.17 393.49 
LS-1.5 0.00 38.69 0.00 38.72 347.12 


The results from Table 8 show that tests SS-13.0-1.0-1.5 and SS-13.0-0.5-1.5 with the same 
sand cover length but different sand cover amounts and thicknesses exhibited similar amount of 
total eroded sediments and runoff yield, with the values being slightly higher for test 
SS-13.0-0.5-1.5 than for test SS-13.0-1.0-1.5 (increments of 72.44 kg and 46.52 m°, respectively). 
However, for tests SS-13.0-0.5-1.5 and SS-6.6-1.0-1.5 with the same amount of sand cover but 
different sand cover patterns, the amount of eroded sediments and runoff yield were largely higher 
for test SS-13.0-0.5-1.5 than for test SS-6.6-1.0-1.5 (increments of 272.73 kg and 108.5 mô, 
respectively). This result indicates that the influences of sand cover pattern on the amount of total 
eroded sediments and runoff yield were much larger than those of sand cover amount. 
Furthermore, the amount of eroded sand was higher for test SS-13.0-1.0-1.5 than for test 
SS-13.0-0.5-1.5 (increment of 65.48 kg), while the amount of eroded loess was higher for test 
SS-13.0-0.5-1.5 than for test SS-13.0-1.0-1.5 (increment of 137.92 kg). This indicates that given 
the same length of sand cover, the thicker sand layer could increase sand erosion, while the 
thinner sand layer could aggravated loess erosion. 

The amounts of eroded loess and eroded sand along with time period during the erosion 
processes in the experimental tests of SS-13.0-1.0-1.5, SS-13.0-0.5-1.5 and SS-6.6-1.0-1.5 are 
shown in Figure 6. The results show large differences in the amount of eroded sand among the 
three experimental tests in the early stage of runoff and sediment generation. And, the amount of 
eroded sand decreased rapidly and tended to become stable during the intermediate and late stages 
of runoff and sediment generation. These results are likely related to the differences in infiltration 
rates between sand and loess layers (Janeau et al., 2003) resulted from discrepancies of relative 
porosities (Van Dijk and Bruijnzeel, 2004). Previous research has shown that the presence of sand 
cover on slopes can extend the time period from rainfall to runoff and sediment generation (4-14 
times) and increase the infiltration rate by 0.2-2.0 times (Xu et al., 2015). The infiltration rate of 
sand layer was greater than that of loess layer, and hence the infiltration was considered as laminar 
flow, which is similar to interflow and may have formed at the sand-loess interface (Eastham et al., 
2000; Hardie et al., 2013). This interflow initiates surface runoff when the water table depth within 
the sand layer reaches a critical value and a peak erosion rate is subsequently attained (lida, 2004; 
Fox and Wilson, 2010). Therefore, the sand layer was transported in the form of a bed load instead 
of a suspended load in the early stage of runoff and sediment generation, and the sand erosion 
showed an obvious peak in the early stage of runoff and sediment generation (Fig. 6b). 


JOURNAL OF ARID LAND 


20 


7 (b) 

Ej = — §S-13.0-0.5-1.5 
z =e — S5-6.6-1.0-1.5 
3 g — §S-13.0-1.0-1.5 
T 

3 3 

le} ao] 

5 3 

Cet v 

9 5 

3 5 

g fo} 

< £ 


1 5 9 13 17 21 25 29 1 5 9 133 17 21 25 29 
Time (min) Time (min) 


Fig. 6 Amounts of eroded loess (a) and eroded sand (b) along with time period during runoff and sediment 
generation in the experimental tests SS-13.0-0.5-1.5, SS-13.0-1.0-1.5 and SS-6.6-1.0-1.5 


By comparing the amounts of eroded loess and eroded sand between experimental tests 
SS-13.0-1.0-1.5 and SS-6.6-1.0-1.5, we found that both the amounts of eroded loess and eroded 
sand were higher for test SS-6.6-1.0-1.5 (one-half slope) than for test SS-13.0-1.0-1.5 (whole 
slope) during the late stage of runoff and sediment generation. Indeed, the time period before 
runoff and sediment generation on SS was much longer than that on LS under the same rainfall 
intensity (Xu et al., 2015). The loess surface of slopes covered by sand with one-half of the total 
length (SS-6.6-1.0-1.5) generated runoff and sediments firstly, but the generation was intercepted 
by the remaining sand surface (one-half slope) and then retained in the sand cover layer, leading 
to a centralized runoff and sediment release after a given time period. This inference was 
confirmed by the runoff coefficients shown in Figure 7. 
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Fig. 7 Runoff coefficients along with time period for experimental tests SS-13.0-0.5-1.5, SS-13.0-1.0-1.5, 
SS-6.6-1.0-1.5, and LS-1.5 


For tests SS-13.0-1.0-1.5 and SS-13.0-0.5-1.5, runoff coefficient being higher than 1.0 occurred 
in the intermediate stage of runoff and sediment generation. For test SS-6.6-1.0-1.5, the runoff 
coefficient increased along with time period and there was an abrupt value (>1.0) in the late stage 
of runoff and sediment generation. Moreover, the runoff coefficient was always lower than 1.0 for 
test LS-1.5. Fox and Wilson (2010) indicated that sand layer coverage changes the generation 
mode of runoff and sediments. Some other studies also demonstrated that the runoff coefficient 
can be greater than 1.0 for sand-covered slopes during runoff and sediment generation 
(Martinez-Murillo et al., 2013; Zhao et al., 2013; Zhang et al., 2017a). With increasing amount of 
sand cover, the time period before runoff and sediment generation increased, demonstrating that 
sand-covered areas at low elevations can intercept runoff and sediments generated from slopes 
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and that sand layers can provide storage and release of runoff and eroded sediments (Ng et al., 
2015). The storage and release processes could be the main reasons why SS generated a large 
amount of eroded sediments. 


5 Conclusions 


The eroded sediments from SS during runoff and sediment generation is dominated by loess, 
which accounts for 83% of the total eroded sediments on average. The sand cover on loess slopes 
transform silt erosion to sand erosion by reducing splash erosion and changing the 
rainfall-induced slope erosion process. The peak sand percentage in eroded sand generally occurs 
in the early stage of runoff and sediment generation, and the proportion of eroded sand from SS in 
the early stage is always higher than that in the late stage. 

Sand covered on loess slopes not only adds additional eroded sediments but also increases the 
amount of eroded loess. The influence of sand cover pattern on runoff yield and the amount of 
eroded sediments is larger than that of sand cover amount. In contrast, given the same length of 
sand cover, a thinner sand layer slightly reduces the sand erosion but significantly aggravates the 
erosion of background sediments. This will lead to a significant increase in the amount of eroded 
sediments from aeolian-fluvial interplay erosion regions to river channels. It should be noted that 
this research mainly focused on the influences of sand cover on erosion processes of loess slopes. 
On the basis of this study, the critical condition of the peak eroded sediments of SS and the 
relationship between erosion energy and eroded sediments can be further studied in the future, 
which are crucial for erosion prediction. 
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